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ABSTRACT 

We study the 850//m emission in X-ray selected active galactic nuclei (AGN) in the ^2 
deg^ COSMOS field using new data from the SCUBA-2 Cosmology Legacy Survey. We find 
nineteen 850/.im bright X-ray AGN in a “high-sensitivity” region covering 0.89 deg^ with 
flux densities of S85o=4-10 mJy. The 19 AGN span the full range in redshift and hard X- 
ray luminosity covered by the sample - 0.7 < z ^ 3.5 and 43.2 < logio(Lx) < 45. We 
report a highly significant stacked 850pm detection of a hard X-ray flux-limited population 
of 699 z > 1 X-ray AGN - S85o=0.71±0.08mJy. We explore trends in the stacked 850pm 
flux densities with redshift, finding no evolution in the average cold dust emission over the 
redshift range probed. For Type 1 AGN, there is no significant correlation between the stacked 
850pm flux and hard X-ray luminosity. However, in Type 2 AGN the stacked submillimeter 
flux is a factor of 2 higher at high luminosities. When averaging over all X-ray luminosities, 
no significant differences are found in the stacked submillimeter fluxes of Type 1 and Type 2 
AGN as well as AGN separated on the basis of X-ray hardness ratios and optical-to-infrared 
colours. However, at logio(L 2 _io / erg s“^ > 44.4), dependences in average submillimeter 
flux on the optical-to-infrared colours become more pronounced. We argue that these high 
luminosity AGN represent a transition from a secular to a merger-driven evolutionary phase 
where the star formation rates and accretion luminosities are more tightly coupled. Stacked 
AGN 850pm fluxes are compared to the stacked fluxes of a mass-matched sample of K- 
band selected non-AGN galaxies. We find that at 10.5 <logio(M*/MQ)< 11.5, the non-AGN 
850pm fluxes are 1.5 — 2x higher than in Type 2 AGN of equivalent mass. We suggest these 
differences are due to the presence of massive dusty, red starburst galaxies in the AT-band 
selected non-AGN sample, which are not present in optically selected catalogues covering a 
smaller area. 
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1 INTRODUCTION 


Studying the connection between star formation in galaxies and ac¬ 
cretion onto the supermassive black-hole is of prime importance 
in un derstanding galaxy formation (e.g. see [Alexander & HickoxI 
I2OI2I for a recent review). The discovery that supermassive black- 
holes are ubiquitous with their mass connected tightly to the stel¬ 
lar bu lge mass of galaxies IMaeorrian et alJl998l : lKormendv & Hoi 
l2013h has led to the realisation that these active galactic nuclei 
(AGN) play a fundamental role in regulating the growth of galaxies. 
Sophisticated galaxy formation simulations now incorporate feed¬ 
back from the AGN, which is shown to be critical in order to ex¬ 
plain t he observed numbe r s of massive galaxies (e .g. lSpringel et al.l 
I2OO5I : ICroton et al.l 1 20061 : IPi Matteo et al.l 1 20051) . In these simu¬ 
lations, massive galaxies are assembled through gas-rich mergers 
which induce luminous bursts of obscured star formation. Star for¬ 
mation and black-hole accretion are fuelled by a common gas sup¬ 
ply and the star formation is initially dust-obscured. As the gas and 
dust are expelled from the galaxy e.g. via AGN-driven winds, the 
accreting black-hole emerges as an X-rayAJV luminous quasar. In 
such a picture of galaxy formation, luminous starburst galaxies, ob¬ 
scured and unobscured AGN are linked through a well-defined evo¬ 
lutionary sequence. Wit hin unified AGN models ( I Antonuccill 1 9931 : 
lUrrv & Padovanilll995h . the obscuration of the central accreting 
source can also be explained purely as a result of orientation effects 
and the distinction between unobscured Type 1 AGN and obscured 
Type 2 AGN is attributed to the line of sight through which they are 
viewed. 


Obscuration in AGN almost certainly results from a com¬ 
bination of orientation and evolution driven factors. One way to 
distinguish between these is to look for evidence for obscured 
AGN having excess star formation in their host galaxies rela¬ 
tive to their unobscured counterparts. Several early studies of X- 
ray luminous AGN using the SCUBA bolometer on the James 
Clerk Maxwell Telescope (JCMT) showed this to be the case 
jPaee et al.ll200lL |2004 Istevens et al.ll2005h . More recently, new 
data at far infrared wavelengths from the Herschel Space Ob¬ 
servatory, has been used to study the connection between star 


Bonfield et alj|201 iklLutz et al.l 20081: Shao et alj 2010l: Page et al. 

20121: iHarrison et al.l 12012 

: iMullanev et all l2012l: iRovilos et al. 

2 OI 2 I: Rosario et alj 20121: 

Santini et al. 20121: Azadi et alj 20141: 


of the 870^m properties of 2 ~ 1 AGN using LABOCA. There 
is now a growing consensus that moderate luminosity AGN tend to 
occupy star-forming host galaxies and that the star formation rate is 
independent of AGN X-ray luminosity in this regime. AGN activ¬ 
ity in galaxies occurs on very different spatial as well as temporal 
scales relative to star formation, which means that the instantaneous 
star formation rate in galaxie s is unlikely to be tig htly coupled to the 
instantaneous accretion rate jHickox et alj2014]) . However, at high 
luminosities, major-merger driven evolution can lead to a tighter 
coupling between AGN activity and star formation. 


Studying the connection between star formation and AGN ac¬ 
tivity at high luminosities requires multi-wavelength data from the 
X-ray to submillimeter wavelengths over reasonably wide fields. 
The COSMOS field is among the most well-studied extragalactic 
fields with a plethora of multi-wavelength photometric and spec¬ 
troscopic data. X-ray data is available from both the XMM-Newton 
and Chandra observatories and the field has been imaged at many 
other wavelengths, in particular with Spitzer and Herschel. More¬ 
over, the wide area covered by many of these datasets now provides 


crucial information for the most luminous, high-redshift galaxies 
and AGN. Recently, the SCUBA-2 bolometer on the JCMT has 
provided a submillimeter map at 850/im over the entir e 2 deg^ field 
as pa rt of the SCUBA-2 Cosmology Legacy Survey l lGeach et al.l 
l2013h . The 850^m data offers several advantages over Herschel 
observations over the same region. Firstly, the superior angular 
resolution and smaller beam size relative to Herschel mean that 
source blending is less of an issue when inferring the average sub¬ 
millimeter properties of AGN (see Fig. [T]l. Secondly, the neg ative 
k-correction at these wavelengths at 2 > 1 l lBlain et alj|2002f) . en¬ 
sures that we are sensitive to roughly constant dust luminosity for 
essentially all high redshift galaxies in an 850/im flux-limited sam¬ 
ple. Furthermore, even in the HerschelSVlRE bands, there can be 
a significant contribution fro m AGN dust heating in the case of the 
most luminous quasars (e.g. lRosario et alj|2012t) . The submillime¬ 
ter dust emission at 850/im is much less prone to such contami¬ 
nation from the AGN. Although it could potentially be more con¬ 
taminated by synchrotron emission, this is unlikely to be an issue 
for the majority of AGN which are radio-quiet The dust emission 
at 850/im can therefore more safely be attributed to heating by star 
formation. At shorter wavelengths it becomes necessary to disen¬ 
tangle the relative contributions of starburst and AGN heating to 
the total infrared flux and linking the two phenomena is therefore 
more complicated, with potential biases arising to derived star for¬ 
mation rates due to different assumptions regarding the intrinsic 
AGN model. 

The aim of this paper is to study the 850/rm properties of 
both unobscured (type 1) and obscured (type 2) AGN in the COS¬ 
MOS field where statistically significant numbers of these popu¬ 
lations now exist. We look for dependences of this submillimeter 
emission on the AGN host galaxy properties and compare to non- 
AGN galaxies. As such, our analysis is intended to provide the first 
SCUBA-2 Cosmology Legacy Survey 850/im view of the X-ray 
AGN population at high redshift. Throughout this paper we assume 
a flat concordance cosmology with Ho=70 km s“^ Mpc“^. 


2 DATA 

2.1 SCUBA-2 850/im 

Our primary dataset used in this analysis is the SCUBA-2 850/im 
map in the COSMOS field obta ined as part of the SCUBA-2 Cos¬ 
mology Legacy Survey (S2CLS: lGeach et alj2013l) . The wide-field 
map covers a total area of ~2 deg^ imaged using a 4x45' PONG 
scanning pattern. The basic principle of the data reduction is to ex¬ 
tract astronomical signal from the timestreams recorded by each 
bolometer in the SCUBA-2 array and map them onto a two dimen¬ 
sional celestial projection. We have used the Dynamical Iterative 
Map-Maker (DIMM) with in the Sub-Millimet re Common User Re- 
duction Facility (SM URF: Ichanin et al.ll20I3l) . We refer readers to 
IChanin et alj l l20I3h for a detailed overview of SMURF, and detailed 
data reduction steps will be given in Geach et al. (in preparation) 
but describe the main steps here. 

Flat-fields are applied to the time-streams using flat scans that 
bracket each observation, and a polynomial baseline fit is sub¬ 
tracted from each bolometer’s time-stream. Each time-stream is 
cleaned for spikes and DC steps are removed and gaps filled. After 
cleaning, the DIMM enters an iterative process that aims to fit the 
data with a model comprising a common mode signal, astronomical 
signal and noise. Next, a filtering step is performed in the Fourier 
domain, which rejects data at frequencies corresponding to angular 
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Figure 1. Stacked 250/im (left), 500/xm (middle) and 850^m (right) images of 2865 radio sources in the COSMOS field from ISchinnerer et alj UoiOl) . The 
superior spatial resolution and angular' scale of the SCUBA-2 data is appai'ent relative to Herschel. For the 850/rm stack, we confirm that the peak emission 
is clearly coming from the central pixels demonstrating that the astrometry of the maps is accurate at the sub-pixel level. The negative “bowling” around the 
bright central source seen in the SCUBA-2 850/rm image, is an artefact of the filtering procedure employed to generate the map. 


scales 9 > 150” and 9 < 2". The next step is to estimate the astro¬ 
nomical signal, which is subtracted from the data. Finally, a noise 
model is estimated for each bolometer by measuring the residual, 
which is then used to weight the data during the mapping process 
in additional steps. The iterative process above runs until conver¬ 
gence is met. In this case, we execute a maximum of 20 iterations, 
or when the map tolerance reaches 0.05. 

The final 850)rm map in the COSMOS field has non-uniform 
depth with typical RMS values of asso ~2 mjy in one half, 
crsso ~4 mJy in the other half and asso ~1 mJy in the central 
region. An 850/im catalogue containing 360 sources has been pro¬ 
duced from this map by searching for all sources above 3.5a in re¬ 
gions with (7850 < 2 mJy (Geach et al. in preparation). This 850/im 
catalogue covers an area of 0.89 deg^. Our search for AGN that 
are individually detected at 850/im covers this smaller area, higher 
sensitivity region. However, when studying the stacked properties 
of the AGN, we utilise the full area in order to get the largest AGN 
samples possible. As stacked fluxes are weighted by the noise (see 
Section S, AGN in lower RMS regions are automatically down¬ 
weighted when calculating average 850/im fluxes. 

We begin by checking the astrometry of this 850/im map by 
stacking a sample of 2865 radio sources in the COSMOS field 
dSchinnerer et al.ir201(]h in this map. The stacking procedure is de¬ 
scribed in more detail in Section|4]and the stacked 850/im image for 
these radio sources is shown in Fig.[T] The peak 850/im emission 
has an offset of 1.4” relative to the radio position. The SCUBA-2 
850/im pixel scale is 2”. This confirms that the maps have the cor¬ 
rect astrometry at the sub-pixel level and can therefore be used in 
our stacking analysis. We also show the stacked images of the same 
radio sources in the Herschel-SPIRE 250 and 500/im bands in Fig. 
[T] The Herschel-SPIRE data in the COSMOS field is taken from the 
Hersc hel Multi-tiered Extragalactic Survey (HerMES: l01iver et al.l 
l2012h . Maps have been produced u sing the Level 2 data p roducts 
from the ESA archive as detailed in ISwinbank et al.l ( l2014h . While 
the stacked 250/im emission appears to be well centred at the ra¬ 
dio position, the 500/rm emission shows an offset relative to the 
radio position. The superior spatial resolution and pixel scale of the 
SCUBA-2 data relative to Herschel-SPIRE, are clearly apparent in 
this plot. The longer wavelength SCUBA-2 850/im observations 
are therefore extremely complementary to already published Her¬ 
schel observations of AGN in this field. 


2.2 AGN Catalogue 

The aim of this paper is to study the submillimeter properties of 
a statistically si gnificant and homo genous sample of AGN. The 
COSMOS field [Scoville et al.l2007h provides the largest and most 
well characterised AGN sample overlapping the new SCUBA-2 
CLS data. The AGN parent sample is a catalogue of XMM-Newton 
detected X-ray point sources, which goes down to a flux limit of 
9.3e-15 erg/s/cm^ in the 2-10 keV band over 90% of the area. 
The flux limit for the faintest sources detected in the catalogue 
are 5e-16 erg/s/cm^ and 2.5e-15 erg/s/cm^ in the 0.5-2 keV and 
2-10 keV bands respectively. We only consider X-ray sources with 
robu st multi-wavelengt h optical and infrared cross-identifications 
from iBrusa et al.l ll201(ll) . which corresponds to ~98% of the X-ray 
sources. We make use of the latest spectroscopic and photomet¬ 
ric re dshifts for these AGN dBongiorno et al.ir2012l : l^lvato et al.l 
I 2 OIII) together with the corresponding spectroscopic and photo¬ 
metric classification!^ Ths Type 1 and Type 2 samples specifi¬ 
cally constitute AGN with broad emission lines in the case of the 
spectroscopic Type Is and narrow emission lines in the case of the 
spectroscopic Type 2s. Eor the AGN with photometric redshifts. 
Type 1 AGN are co nsidered to be those fit by templates 19-30 in 
ISalvato et al.l ( 1201 ih whereas Type 2 AGN are those fit by tem¬ 
plates 1-18 or >100 (which correspond to galaxy templates). The 
photometric redshift accuracy is in general very good aAz/(I + z) 
= 0.014 for Iab < 22.5 and 0.015 for Iab < 24.5 dSalvato et al.l 
l2009h . Photometry at mid infra-red wave lengths is from the Spitzer 
S-COSMOS Survey dSanders et al.l2007h . There are a total of 1797 
AGN with multi-wavelength cross-identifications in these cata¬ 
logues. The observed hard-band (2-10 keV) X-ray luminosity ver¬ 
sus redshift for both the Type 1 and Type 2 AGN can be seen in 
Fig. 1^ These X-ray l uminosities have b een calculated from the 
hard-band fluxes from iBrusa et^ dTOlOh assuming a photon in¬ 
dex of r=1.8. In the case of the Type 1 AGN that are detected in 
the soft (0. 5-2 keV) band but not in the hard band, the soft X-ray 
fluxes from iBrusa et al.l d201(ll) are used to estimate the hard X-ray 
flux, once again assuming r=1.8. We note that no absorption cor¬ 
rection has been applied to these luminosities as the X-ray data 
for the entire sample is in general of insufficient quality to allow 
full spectral fitting. The Type 1 AGN have negligible absorption 


^ http://www2011 .mpe.mpg.de/XMMCosmos/xmm53_release/ 
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and in the case of the Type 2 COSMOS AGN that are detected in 
the hard X-rav. lLusso et al ] |20ll|) estimate the average shift in the 
hard X-ray luminosity introduced by an absorption correction to be 
< Alog(L 2 -io) >= 0.04 ± 0.01. In Appendix A we compare 
the hard X-ray lum i nositi es derived in this paper to those derived 
bv iBrightman et alj ( l2014h for a subset of the COSMOS AGN with 
Cluindm data and where complex absorption spectral models can 
be fit, demonstrating that there are no significant biases in the X-ray 
luminosities used in this paper. 

In the majority of the analysis presented in this paper we 
consider the average stacked submillimeter fluxes of AGN as 
a function of various properties. For these stacking studies, we 
construct flux-limited populations of AGN restricted to sources 
with S>2.5xl0“^® erg s“^ cm~^ in the hard (2-10 keV) band 
JLusso et alj201lh . As previously mentioned. Type 1 AGN that are 
detected in the soft band but not in the hard band, have had their 
hard band fluxes estimated from the soft X-ray flux using a T = 1.8 
power-law SED. Type 2 AGN that are detected in the soft band only 
have been discarded for the majority of the stacking analysis due 
to their essentially unknown dust columns and therefore intrinsic 
luminosities. These soft band detected Type 2 AGN are however 
included when searching for AGN that are individually detected at 
850/rm in Section[3 The hard-band X-ray flux limit is illustrated in 
Fig.l^and corresponds to Lx ^ lO'^^ erg/s at 2 > 1, where the X- 
ray luminosity should be completely dominated by AGN activity. 
We also restrict ourselves to only those AGN at 2 > 1 for all the 
stacking analysis. At 2 > 1, the 850/im flux is essentially invariant 
with redshift for a given dust luminosity and therefore star forma¬ 
tion rate, allowing us to average the submillimeter fluxes of AGN 
over relatively wide redshift bins in order to improve the statistical 
significance of our results. The 850/rm data does not however trace 
the peak of the dust SED at these redshifts and there could poten¬ 
tially be biases in how effectively this submillimeter flux is tracing 
star formation due to variations in the template SEDs over this rest- 
frame wavelength range. In Section l^TI we will therefore explicitly 
check the dependence of 850^m flux on star formation rates. 

Finally, there is a concern that the 850/im flux could be con¬ 
taminated by synchrotron emission for populations of radio-loud 
AGN. The fraction of radio loud AGN in X-ray selected AGN sam¬ 
ples i s expected to be small (<5% for Type 1 AGN from iHao et al.l 
l20Ti) . Nevertheless, we match the X-ray AGN sample to the VLA 
radio catalogue dSchinnerer et alj [201011 and remove all sources 
with Si.4GHz >0.5 mJy. This removes 10 Type 1 and 7 Type 2 
AGN at 2 > 1 from the sam^e, which represents ~2% of both 
the Type 1 and Type 2 sampleo The radio flux limit imposed cor¬ 
responds to a synchrotron 850/rm flux of <0.01 mJy assuming a 
synchrotron spectral index of 0.7. Therefore we can be confident 
that the 850/im stacks that we study are completely dominated by 
thermal dust emission. Our final z > 1 AGN sample used in the 
stacking analysis, totals 699 AGN. These are comprised of 428 
Type 1 AGN (314 spectroscopic, 114 photometric) and 271 Type 
2 A GN (46 spectroscopic, 2 25 photometric). 

iBongiomo et al.l ( 1201 2h have used spectral energy distribution 
(SED) fits to the multi wavelength COSMOS photometry, in or¬ 
der to estimate stellar masses, star formation rates and other host 
galaxy parameters for these Type 1 and Type 2 AGN. Both galaxy 
and AGN components are fit allowing us to study the dependence 


^ Including these radio-detected AGN in the 850/rm stacking does not 
change the average 850^m flux derived for our sample in Section|4]to within 
0.1 mJy 


of 850^m flux densities on a wide range of AGN and host galaxy 
properties. These SED fits incorporate photometry at optical, near 
infrared as well as mid infrared wavelengths. Around 80 per cent 
of the sources are detected in the Spitzer MIPS 24/im band. In ad¬ 
dition, several of the AGN are also detected in the MIPS IQpm 
band as well as in the Herschel-PACS lOO^m and 160/rm bands 
dBongiomo et al.ir2012h . All of this multi-wavelength photometry 
is used in the SED fitting and both pure AGN and pure galaxy tem- 
pl ates are allowed in the fits. The galaxy templates constitute a set 
of iBruzual & Charlod (l2003h stellar population synthesis models 
with a Chabrier ( l2003h initial mass function and assuming expo¬ 
nentially declining star formation histories with varying amounts 
of d ust reddening. The A GN template is the mean quasar SED 
from iRichards et all (l2006h . again with varying amounts of redden¬ 
ing. This choice of AGN template is appropriate for the luminous 
XMM -COSMOS AGN investigated in this work. IBongiomo et al.l 
( l2012h have already discussed the accuracy of the stellar mass esti¬ 
mates for both the Type 1 and Type 2 AGN samples and when con¬ 
sidering the properties of the AGN as a function of stellar mass, we 
only consider those AGN that have robust stellar masses as spec¬ 
ified by the flags parameter. As expected, > 99.5 per cent of the 
Type 2 AGN have robust stellar mass estimates. 


2.2.1 Limitations of the X-ray AGN Catalogue 

In this analysis we consider the submillimeter properties of a hard 
X-ray flux limited catalogue of AGN. However, before proceed¬ 
ing it is important to highlight the selection biase s and incom¬ 
pleten esses of an AGN catalogue selected in this wav. lDonlev et al.l 
( l2012h have used a mid infrared IRAC selection for AGN candi¬ 
dates in the COSMOS field demonstrating that only ~38% of these 
AGN candidates have X-ray counterparts due to the relatively shal¬ 
low flux limit of the XMM-COSMOS hard X-ray data. Hence, we 
may potentially be missing a significant population of heavily ob- 
scured to mildly Compton thick AGN in this analysis. Recently 
iLacv et alj ( 12015h have also used mid infra-red selected AGN sam¬ 
ples to demonstrate that X-ray surveys begin to become incomplete 
at 2 > 1.6 where a deficit in the X-ray luminosity function is seen 
relative to the mid infra-red luminosity function even in the highest 
luminosity bins. We refer readers to these papers for a detailed con¬ 
sideration of the incompleteness in hard X-ray selected samples of 
AGN. 


2.3 Non-AGN Galaxy Catalogue 

Finally, as a control sample to compare the AGN properties to, we 
also util ise a galaxy sample from the Ultra VISTA survey Data Re¬ 
lease 1 dMcCracken et al.ll2012h over the s ame region. We wor k 
with the TC-band flux limited sample from iMuzzin et all ( l2013ll . 
which essentially represents a mass-limited sample of galaxies. 
Corresponding photometric redshifts, stell ar masses and othe r 
SED-fitting parameters are also taken from iMuzzin et al.l (l2013l) . 
An advantage of using a galaxy catalogue that is selected in the 
redder Tf-band, is that these catalogues are potentially more sensi¬ 
tive to dusty galaxies that would be detected at 850/im compared 
to galaxy catalogues constructed using data at optical wavelengths. 
We note that several previous studies that we compare our results 
to in this paper, have used non-AGN galaxy samples selected using 
flux-limited surveys at bluer (optical) wavelengths. 
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Figure 2. Redshift versus observed hard (2-10 keV) band luminosity for Type 1 AGN (left) and Type 2 AGN (right). Both spectroscopic and photometric 
AGN are shown. The circles mark actual detections whereas the downward triangles represent upper limits in the case of the Type 2 AGN and hard-hand 
luminosities estimated from the soft-band flux in the case of the Type 1 AGN. The solid lines denote the flux limit applied to the data to generate homogenous 
populations of the two populations for our stacking analysis. The vertical dashed lines correspond to z > 1 which is the redshift cut applied before stacking 
such that the 850^m flux is invariant with dust luminosity as a function of redshift. 


3 850/rm DETECTED AGN 

We begin by looking for AGN that are detected at >3.5a in the 
SCUBA-2 CLS 850/im catalogues. As stated above, only regions 
with (T 850 < 2mJy were used in the construction of the 850/rm 
catalogue, which therefore does not cover the entire field. A match¬ 
ing radius of 7.5" is used to associate 850^m catalogued sources 
with our AGN sample. The AGN positions in a ll cases correspond 
to the optical positions from iBrusa et al.l jlOlOl) . which have much 
smaller uncertainties than Ihe XMM-Newton X-ray positions. Given 
the errors on the optical positions are likely to be negligible, we fol¬ 
low |ivison|er^ (( 2003 ) and assume that for a signal-to-noise ratio 
(SNR) of 3.5, the positional uncertainty, (jpos=0.6xFWHM/SNR. 
For a FWHM of 14.5" at 850/rm, o-pos is therefore 2.5" and the 
7.5" matching radius corresponds to 3apos and is therefore ex¬ 
pected to encompass >99.7% of true counterparts. 

There are a total of 19 850pm sources that lie within 7.5" of 
an X-ray AGN and there are no instances of multiple submillimeter 
sources matched to the same X-ray AGN. All 19 matches are pre¬ 
sented in Table[T]and we estimate the corrected Poisson probability, 
Pcorr of these being chance associations based on the background 
density of all X-ray AGN that have fluxes brighter than the target 
being considered. Specifically: 

Praw = 1 - exp (•7rr^N(> S)) 

Pdet = 1 - exp (7rri^^N(> Slim)) 

Pcorr = 1 - exp ^Praw + In ^ ^ 

where r is the separation between the submillimeter source and 
the X-ray AGN, rmax is the maximum search radius and N(>S) 
and N(>Siiin) represent the background surface density of AGN 
brighter than the source being considered and brighter than the 
flux limit of the survey respectively. We estimate N(>S) us- 
ing the source c ounts in the soft and hard X-ray bands from 
ICappelluti et alj l l2009l) and assume Siim=5e-16 erg/s/cm^ and 
2.5e-15 erg/s/cm^ in the 0.5-2 keV and 2-10 keV bands respec¬ 
tively. Given the low surface density of X-ray AGN and submil¬ 
limeter sources, all 19 matches have very small probabilities (< 


1%) of being chance associations (see Table [T}. The 19 sources 
are comprised of 6 spectroscopically confirmed broad-line AGN, 1 
spectroscopically confirmed narrow-line AGN, 3 AGN with pho¬ 
tometric redshifts that are best fit by Type 1 AGN templates and 
9 AGN with photometric redshifts that are best fit by Type 2 
AGN or galaxy templates. Thus out of a sample of 360 850pm 
bright sources, we find that -^5% are associated with X-ray lu¬ 
minous AGN. The XMAf-COSMOS data used in this work cor¬ 
responds to re latively bright AG N X-ray flux limits as can be 
seen in Fig. [2l IWang et alT ( 1201 3h identify the X-ray counterparts 
to ALMA detected submillimeter bright galaxies in the Extended 
Chandra Deep Field South (E-CDFS). Chandra observations are 
used to go down to X-ray luminosities of ~7x 10"^^ erg/s and the 
authors find a larger AGN fraction of ~17% among the submil¬ 
limeter galaxies. However as can be seen in their figure 10, this 
AGN fraction decreases as a function of inc reasing X-ray lumi ¬ 
nosity. Down to similar X-ray luminosities as IWang et ^ (l2013l) . 
ISvmeonidis et al.l (1201 4l) find an AGN fraction of 18% among in- 
frar ed luminous galaxies in the CDF-N and CDF-S fields at z < 
1.5. [Alexander et ^ ( l2005h have used ultra-deep X-ray observa¬ 
tions in the CDE-N to study the AGN fraction in Ssso > 4 mjy 
sub-mm bright galaxies at z > 1, finding that the majority host 
AGN activity. However, the X-ray data from Clumdra used in that 
work, is an order of magnitude fainter than our XMM-COSMOS 
sample. In Fig [3 we plot the redshift versus hard (2-10 keV) X- 
ray lumin osity for the submillim eter detected AGN in this work 
as well as [Alexander et al.1 ( l2005t) an d lWang et al.l ( l2013h . We also 
show here the subm illimeter detected X-ray absorbed quasars from 
IStevens et al.l ( l2005h . which were selected using X-ray surveys with 
ROSAT, XMM-Newlon and Chandra and followed up at sub mil- 
limeter wavelengths using SCUBA. The IWang et ^ ( 1201 3l) ob¬ 
served 0.5-8 keV luminosities have been converted to 2-10 keV 
luminosities assuming the intrinsic photon indices derived for each 
o bject in that p a per w hereas conversions from 0.5-2 keV fluxes 
in IStevens et al.l ( 120051) assume F = 1.8. Our work can immedi- 
ately be place in the co ntex t of these previous studies. While both 
[Alexander et alj ( [2005[) and [Wang et alj ( [20131) probe significantly 
fainter in terms of X-ray luminosity, our sample fills in the luminos¬ 
ity gap between these studies and the individual follow-up obser- 
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Figure 3. Redshift versus observed hard X-ray (2-10 keV) luminosity for 
the X-ray AGN that are individually detected at submillimeter wavelengths 
in this work and compared to individua lly d etected submillime ter bright 
X-ray AGN from [Alexander et ^^_j200^_andj^^nge^^ ^2013^ and sub¬ 
millimeter detected quasars from IStevens et al. ilOOSh . Downward triangles 
denote upper limits. The hard X-ray flux limit and z > 1 redshift limit ap- 
plied prior to our stack ing a nalysis, are denote d by the dashed lines. Both 
[Alexander et al.l i2005l) and IWang et^ j2013h probe fainte r X-ray lumi¬ 
nositie s compared to our work, while the X-ray AGN from IStevens et all 
i2005l) that were individually followed up, probe very high X-ray luminosi¬ 
ties. Our study is useful in bridging the gap between these different datasets. 


vations of X-ray bright quasars bv IStevens et al.l ( l2005h . The X-ray 
brightest sub-mm AGN in our sample is brighter than any o f the 
X-ray AGN detect ed in the suh-mm by lAlexander et all ( l2005h and 
IWang et ^ l l2013h . It corresponds to XID18 and is an extremely 
red A GN with evidence for outflowing gas with velocities of ~300 
km/s terusa et al.ll2015h . This class of objects will be discussed 
later in Section IdAl 

Herschel fluxes f or all SSOttm de t ected sources are extracted 
from the catalogues of ISwinbank et al.l l l2014h by matching the op¬ 
tical positions of the AGN in our catalogue to the 24/rm positions in 
the Herschel catalogues, using a matching radius of 3". The Her¬ 
schel catalogues have alrea dy been corrected for blending using 
the 24^m positional priors Jswinhank et al.ll2014h . These fluxes 
are presented in Appendix B together with other derived properties 
for these submillimete r bright X-ray AGN. A s the Herschel fluxes 
haye been extracted by lSwinbank et al] l l2014h using the 24/im po¬ 
sitional priors, all 24^m detected sources haye measured Herschel 
fluxes whereas sources not detected at 24^m are not in the Herschel 
catalogues. 






































Table 1. Summary of X-ray AGN with potential 850/im bright identifications in the COSMOS field. The con'ected Poisson probability of a chance association 
is quoted as calculated from the soft X-ray flux and the hard X-ray flux respectively. Sources with photometric redshifts are mai'ked with a pz. Sources with 
hard X-ray luminosities marked with an asterisk are Type 1 AGN that are only detected in the soft band and where the hard band luminosity has been calculated 
assuming P = 1.8. 


XID 

RA 

Dec 

Redshift 

loglo(L2-lo) 
erg s“^ 

Separation 

// 

So.5-2keV 
erg s“^ cm“^ 

P0.5-2keV 

S2-10keV 
erg s“^ cm“^ 

P2-10keV 

Ssso 

mjy 

13 

150.00924 

2.2755119 

0.850 

43.95 

3.42 

2.11e-14 

0.0006 

2.9e-14 

0.0009 

4.5±1.1 

18 

150.13304 

2.3032850 

1.598 

44.92 

5.92 

1.80e-14 

0.002 

6.04e-14 

0.0009 

4.4±1.2 

139P^ 

150.04184 

2.6294839 

0.739 

43.69 

3.30 

9.52e-15 

0.001 

2.22e-14 

0.001 

7.4±1.7 

160 

150.15839 

2.1396031 

1.825 

44.28 

1.48 

2.39e-15 

0.001 

1.02e-14 

0.0008 

8.1±1.2 

246P^ 

150.05100 

2.4938550 

2.342 

44.56 

3.65 

8.81e-16 

0.008 

1.09e-14 

0.003 

6.4±1.7 

250 

150.06455 

2.3290511 

2.446 

44.14* 

3.04 

2.41e-15 

0.004 

<3.78e-15 

- 

4.0±1.1 

270P^ 

150.10612 

2.0144781 

1.883 

<44.13 

1.83 

1.17e-15 

0.003 

<6.64e-15 

- 

9.8±1.7 

278P^ 

150.09308 

2.1014100 

2.269 

<44.46 

2.61 

1.46e-15 

0.004 

<9.39e-15 

- 

5.0±1.4 

353 

150.08438 

2.2904881 

1.112 

43.70 

2.04 

2.34e-15 

0.002 

8.57e-15 

0.002 

4.7±1.0 

402P"^ 

150.25226 

2.2619131 

1.078 

43.57 

2.86 

2.22e-15 

0.003 

6.82e-15 

0.003 

5.7±1.5 

415P"’ 

149.96981 

2.1834869 

1.558 

43.33* 

2.90 

1.04e-15 

0.005 

<8.31e-15 

- 

5.5±1.4 

469P2 

150.09685 

2.0215000 

3.362 

<44.45 

2.22 

1.08e-15 

0.004 

<3.73e-15 

- 

6.6±1.7 

10675 

150.19262 

2.2198489 

3.090 

43.92* 

0.28 

8.45e-16 

0.0001 

<1.73e-15 

- 

5.7±1.2 

10809P"' 

150.20758 

2.3816111 

1.288 

<44.15 

1.79 

1.12e-15 

0.003 

<1.7e-14 

- 

6.5±1.5 

30182P^ 

150.12541 

2.6978931 

0.742 

43.23 

6.00 

6.26e-15 

0.005 

7.52e-15 

0.009 

6.9±1.7 

53922 

150.09454 

2.7029131 

0.850 

43.38 

1.81 

2.25e-15 

0.002 

7.74e-15 

0.001 

8.6±1.8 

54440PZ 

150.06308 

1.9446778 

1.947 

43.91 

4.39 

1.13e-15 

0.009 

<3.74e-15 

- 

7.7±1.7 

60070P^ 

150.22719 

2.2324781 

2.066 

<44.01 

1.93 

7.33e-16 

0.003 

<4.11e-15 

- 

5.8±1.4 

60490^"^ 

150.10544 

2.1852681 

1.089 

43.42 

5.33 

<7.57e-16 

- 

4.76e-15 

0.011 

6.9±1.1 




4 AVERAGE 850^lm EMISSION EROM STACKED 
SUB-SAMPLES 

The majority of the X-ray AGN are not individually detected in the 
850^m catalogue, which goes down to a 3.5<j flux limit of ~3.5 
mjy in the highest sensitivity regions. To understand the average 
submillimeter properties of the AGN we therefore have to study 
stacked 850^m maps. We begin by constructing an inverse vari¬ 
ance weighted stacked image at 850/im for all 699 z > 1 AGN 
(428 Type 1 and 271 Type 2) that will be used from hereon in our 
stacking analysis. Specifically: 


^ X /i ;./(4 X 4 ) 

where Sij represents the stacked 850/im image, Pij is the point re¬ 
sponse function of SCUBA-2 at 850/im generated by stacking all 
>10(7 sources in the map, fij are the 40 x 40 pixel (80x80") flux 
images of each of the N sources going into the stacks and aij are 
the corresponding RMS images. Stacked images generated in this 
way can be seen in Fig. |4] for both the AGN sample and a set of 
random points in the map, where any random points within 7.5" 
of a catalogued 850/im source, have been removed. Fig. |4] clearly 
demonstrates the presence of a statistically significant 850/im sig¬ 
nal from the z > 1 X-ray AGN relative to the random points. The 
typical RMS noise in the maps is well above the confusion limit. 
As stated in Section [ZT] the noise is non-uniform over the full 2 
deg^ map. Out of the 699 z > 1 sources being stacked, 80% lie 
in regions with RM S noise <3 . 5 mJy . From the 850/im source 
counts presented in ICasev et^ (l2013h . we estimate that the in¬ 
tegral source count at this RMS noise level is only ~630 deg~^. 
Assuming a beam FWHM of 14.5" for SCUBA-2 at 850/im, the 
beam density is ~78,000 deg“^. Hence only 1 in 120 beams are 
expected to be affected by confusion at this RMS noise. In the low¬ 
est RMS regions in our map ((j=lmJy), the integral source counts 
reach ~3000 deg“^ and 1 in 25 beams are affected by confusion, 
which is still a small effect. 

Before investigating the dependence of the 850/im emission 
on the AGN properties, we first describe our methodology for cal¬ 
culating the stacked 850/im fluxes. Throughout this analysis the 
stacked 850/im flux is measured as the median or mean value of 
the central pixel in 100 bootstrapped inverse variance weighted 
stacked images (see Eq. [2]l. We have already checked in Section 
12. II that the astrometry of the SCUBA-2 map is accurate to within 
a pixel so this central pixel flux provides the most conservative and 
unbiased estimate of the average flux. The central pixel flux is mea¬ 
sured to be 0.71 ±0.08 mJy for all 699 AGN. The median redshift 
of the sample is 2 = 1.66 and the median hard X-ray luminosity 
is logio(L 2 -io) = 44.2. We emphasise that these X-ray luminosi¬ 
ties correspond to the most luminous AGN in th e “quasar” regime . 
For comparison, the average 870/jm flux of th e iLutz et aH ( 1201 Oh 
AGN is 0.49±0.04mJy. The lUutz et al] ( l2010h sample has both a 
lower median redshift (z = 1.17) and a lower median hard X-ray 
luminosity - logio(L 2 -io) = 43.6 calculated from the median hard 
X-ray flux in that paper and using our assumed value of F = 1.8 
as described in Section Given these differences in the median 
properties of the two samples as well as the correction from 870/im 
to 850/im for a typical single temperature greybody SED and the 
statistical errors in b oth measurements , our results are broadly con¬ 
sistent with those of iLutz et al.l (l201(lh . Excluding the individually 
detected 850/im X-ray AGN described in Section[3from the stacks. 
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Figure 4. 80 x 80” weighted, stacked 850/im images centred on the posi¬ 
tions of the 699 2 > 1 X-ray AGN in our flux limited sample (top) and 
centred on exactly the same number of random points (bottom). The AGN 
stack includes 428 Type 1 AGN and 271 Type 2 AGN. The colour scaling 
is the same in both images. The contours are at l-5cr and the solid contours 
correspond to positive fluxes whereas the dashed contours correspond to 
negative fluxes. We clearly observe a statistically significant 850/tm detec¬ 
tion from the AGN relative to the random points. 


we find that the average flux falls to 0.60±0.07 mJy, which is still 
a >8 (t detection. 

We have checked that the median and mean stacked fluxes are 
very similar in both cases. The error on the flux is calculated using 
100 bootstrap samples. We also cal culate stacked flu xes using the 
publicly available code SI MS TACK dViero et alj|2013l) . The code is 
designed to correctly account for clustering of sources within the 
scale of the SCUBA-2 beam which is assumed to be 14.5". The 
method essentially involves performing a regression of the true flux 
map with a ‘hits’ map, where the ‘hits’ map is created by counting 
the number of sources to be stacked that contribute to each pixel in 
the map. We stack the 699 flux-limited 2 > 1 AGN with the 2 < 1 
AGN. If the photometric redshifts are accurate, there should be no 


































clustering between the high and low redshift AGN populations. The 
stacked 850/rm flux for the 699 z > 1 AGN is 0.82±0.02 mjy from 
SIMSTACK which is slightly higher than our central pixel flux mea¬ 
surement. The higher SIMSTACK flux is likely due to the fact that 
the peak 850^m emission in Fig.|4]is offset from the central pixel. 
The SIMSTACK errors reflect the regression errors and do not ac¬ 
count for any errors due to sample variance. Our bootstrap errors 
are therefore more conservative. The SIMSTACK results demon¬ 
strate that the central pixel flux represents a conservative estimate 
and is not biased up due to clustering. For the rest of the analy¬ 
sis, we therefore use the central pixel flux to measure the average 
submillimeter emission for all the sub-samples we construct. As 
a final check, given that some regions of the 850pim wide S2CLS 
COSMOS maps are significantly noisier than others (Section fT. 11 1 , 
we also restrict our stacks only to the regions with crsso < 2 mJy to 
mitigate against the effects of flux boosting. There are 341 z > 1 X- 
ray AGN in our flux-limited sample in these high sensitivity regions 
and the mean flux is S85o=0.69±0. 10 mJy, consistent with the re¬ 
sults derived using the entire 850/rm wide map. We will therefore 
proceed with using the 850/rm data over the full COSMOS wide 
area for all the stacking analysis that follows. 

Before continuing, it is important to highlight the assumptions 
inherent in our work. Throughout this paper, we assume that the 
850^m flux is tracing the total amount of cool dust emission in 
the AGN host galaxies and that this cool dust is primarily heated 
by star formation and has a negligible contribution from the AGN. 
Given that the 850/rm flux is invariant with redshift at 2 > 1, we 
can average these fluxes over relatively broad redshift bins without 
bias. Note that this is not the case for Herschel fluxes for exam¬ 
ple, which trace the peak of the dust SED at 2 ~ 1 — 3, and are 
expected to be varying more strongly with redshift. As we do not 
have large enough AGN samples to be able to split our AGN into 
narrow redshift bins, we restrict our analysis to the 850/^m prop¬ 
erties only rather than constructing full SEDs from the Herschel 
and 850/rm data. The 850^m emission can also only be used as 
a proxy for the star formation rate provided the dust temperature 
is not varying significantly between AGN sub-samples, and also 
assuming that the dust heating at these wavelengths is dominated 
by recent star formation rather than being dominated by cold dust 
heate d by old stars and the interstellar medium (e.g. iBourne et al.l 
l2013h . We therefore begin by testing how the 850/im flux varies 
with star formation rate. 

4.1 Star Formation Rate 

iBongiomo et alj ( 1201 2l) have estimated the star formation rates 
for the COSMOS AGN using SED-fitting to the available multi¬ 
wavelength photometry dBrusa et al.ll201^ - SERsed hereafter. 
These can be used to check the dependence of the 850^m flux den¬ 
sity on star formation rate. As described in Section |2^ the SED 
based star-formation rates include photometry in the Spitzer IRAC 
as well as MIPS 24/rm bands in most cases. In some cases, the SED 
fits also include longer wavelength photometry at MIPS 70/im as 
well as at Herschel 100 and 160^m. We select only the Type 2 
2 > 1 AGN with robust star formation rate es t imate s as speci¬ 
fied by the flags parameters in iBongiorno et al.l (l2012h . We note 
that soft X-ray detected Type 2 AGN that are not detected in the 
hard band are also included in this analysis to improve the statis¬ 
tics. Here we are primarily concerned with checking whether the 
850^m fluxes correlate with other independently derived star for¬ 
mation rate constraints and knowledge of the true dust columns and 
therefore intrinsic X-ray luminosities of the AGN is not relevant. 


In luminous, unobscured Type 1 AGN, the big blue bump typically 
seen in quasar spectra may mimic the UV bump from massive stars 
so star formation rates based on SED fitting are no longer reliable. 
Hence Type 1 AGN are not considered in this test. 

We split the Type 2 AGN into two SFRsbd bins at 
SFR>lOM 0 yr“^ and SER<lOM 0 yr“^. The median redshifts are 
2 = 1.36 and 2 = 1.77 for the low and high SERsed sub-samples 
respectively. The median star formation rates are 1 . 5 M 0 yr~^ and 
33 M 0 yr“^ while the mean star formation rates are l.lM 0 yr“^ 
and 46M0yr“^ in the two bins. The average 850^m fluxes for 
these sub-samples are 0.3±0.1 miy and 1.0±0.2 mJy respectively 
and are therefore clearly higher in the higher SFR bin. The ex¬ 
cess 850/rm flux in the high SFR bin is significant at the ~3 (j 
level. Several previous studies have instead used the 250/rm fluxes 
as a proxy for SFR. We therefore also calculate average Herschel 
250^m fluxes for the Type 2 AGN in these two SFR bins although 
we caution that the bins encompass AGN over a relatively broad 
range of redshifts, which could be adding scatter to the 250^m 
fluxes. The 250/im fluxes are calculated by injecting the AGN in 
each SFR bin in turn into the Herschel SPIRE 250/rm map from 
ISwinbank et al.l ( l20I4h . which has had all 24/im detected sources 
removed. The maps are then stacked at the AGN position and 
the average flux is once again measured as the value at the cen¬ 
tral pixel. Note that we have also checked the astrometry of the 
250^m map in Section ITT] The 250^m fluxes in these two SER 
bins are 2.0±0.3 mIy and 5.3±0.6 mIy respectively. In order to 
check whether the 250/im and 850^m stacks are consistent with 
each other, we fit a single temperature greybody with /3=2.0 to 
these data points in the two SFR bins. The dust temperatures are 
24±6K and 27±5K respectively, consistent with a single temper¬ 
ature in both SFR bins. From the integral under the greybody, we 
can then calculate the far infrared luminosity (between 60-300/rm 
to eliminate any contributions from the AGN) and average star for¬ 
mation rate for each stack. These are logio(LFiR/L 0 )=lO. 9 ±O .3 
corresponding to SFR= 12 ± 6 ^M 0 yr“^ for the low SFR stack and 
logio(LFiR/L 0 )= 11 . 5 ±O .2 corresponding to SFR= 6 O± 2 oM 0 yr“^ 
for the high SFR stack. In both cases the star formation rates de¬ 
rived from the Herschel a nd 850/rm data are high er than the mean 
star formation rates from IBongiorno et al.l l l20I2h although the er¬ 
ror bars and scatter in both average SFR measurements are large. 
In Appendix B where we have shown best-fit far infrared SEDs 
for our X-ray AGN that are individually detected at 850pim, we 
also show that the far infrared luminosity and star formation rate 
as measured from these best-fit SEDs, correlates well with both the 
250 and 850/im fluxes (Eig. IB2t . In Appendix B, the dust mass 
from full SED fitting has also been plotted against both the 250 and 
850^m fluxes and in general shows a weaker correlation with the 
infra-red and sub-mm flux densities compared to the far infra-red 
luminosity. We will therefore proceed by assuming that the stacked 
850/rm flux in our z > 1 AGN is giving us a direct estimate of the 
average star formation rate, with the caveat that for galaxies with 
low levels of star formation, there may also be a contribution to the 
dust heating at 850/im from old stars and the interstellar medium. 

4.2 Redshift 

The size of our current sample necessitates stacking over relatively 
broad redshift bins. Evolutionary effects could therefore influence 
any derived results and introduce differences in stacked 850/rm 
fluxes between different sub-samples. Before looking for trends 
with other AGN properties we therefore first consider the redshift 
evolution of the 850^m fluxes over the full redshift range consid- 
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Figure 5. Average 850^m flux versus redshift for both Type 1 and Type 2 Figure 6. Average 850/rm flux versus hard X-ray luminosity for both Type 
AGN. 1 and Type 2 AGN. 


ered in our stacked analyses. In order to do this, we split the z > 1 
Type 1 and Type 2 AGN into three redshift bins and consider their 
stacked 850/rm fluxes. The results are illustrated in Fig. Both 
the Type 1 and Type 2 AGN submillimeter fluxes show little evolu¬ 
tion with redshift over the redshift range probed in this work. In all 
three redshift bins, the Type 1 and Type 2 AGN 850^m fluxes are 
consistent with each other given the error bars. When considering 
trends in AGN submillimeter fluxes with luminosity, obscuration 
and AGN activity in the following sections, we will therefore as¬ 
sume that the 850^m fluxes do not evolve with redshift over the 
redshift range considered. 


4.3 X-ray Luminosity 


We now explore trends in these 850pim fluxes with the hard X-ray 
luminosity, which can be used as a proxy for the energy emitted 
by the central accreting power source. Several recent studies have 
found a reasonably flat r elationship between AGN luminosity and 
star formation rate (e.g. Hamson^t al.l l2012l : iRosario et alj|201^ : 
lAzadi et^l2014l : IStanlev et al.ll2015h . At modest X-ray luminosi¬ 


ties where large samples of X-ray AGN with far infrared and sub¬ 
millimeter data now exist, these results are now robust. However, 
the situation is le ss clear at very high X-ray luminosities. While 
IPage et al] ( 1201 2h have suggested that the far infra-red fluxes and 
therefore star formation rates of AGN with high X-ray luminosi¬ 
ties, ar e su ppressed relative to A GN of lower luminositv. lLutz et al.l 
(l2010h and iRosario et al.l (1201 2h instead And an increase in star for¬ 
mation rate at high X-ray luminositi es, with the results be ing more 
pronounced for low redshift AGN in iRosario et al.l ( l2012h . The re¬ 
sults from recent hydrodynamical simulations suggest that there is a 
large scatter in the Lx-SFR relation ove r the dynamic range probed 
by current surveys dSiiacki et al.ll2014h . which may explain some 
of the discrepancies between previous studies. 

The SCUBA-2 photometry presented in this paper now al¬ 
lows an independent investigation of the relation between star for¬ 
mation and AGN luminosity with the advantage that the submil¬ 
limeter flux is less contaminated by AGN emission than the Her- 
schel fluxes in the case of the most luminous quasars. As illus¬ 
trated in Fig. 12 at 2 > 1, the COSMOS field is only complete at 
relatively high X-ray luminosities. We select Type 1 and Type 2 


AGN in three X-ray luminosity bins - 43.5< log]^Q(L 2 -io) < 44, 
44.0< logjQ(L 2 -io) < 44.4 and log j^q(L 2 -io) > 44.4. The AGN 
that are individually detected at 850/rm (see Section [2 are also 
included in the stack. There are only 7 X-ray AGN out of the 19 
that satisfy both the hard X-ray flux limit and the z > 1 redshift 
cut and we have checked that inclusion of these AGN does not 
change the stacked fluxes quoted within the error bars. Note that 
the lowest X-ray luminosity bin is incomplete at 2 > 2 so the re¬ 
sults here should be interpreted with caution. The average 850^m 
fluxes are presented in Tableland illustrated in Fig. [6] While the 
Type 1 submillimeter flux appears to be fairly independent of X-ray 
luminosity, there is some evidence that the most X-ray luminous 
Type 2 AGN have higher submillimeter fluxes. The Type 2 AGN at 
logj^Q(L 2 -io) > 44.4 have 850/rm fluxes that are ^2a higher than 
in the previous X-ray bin although the numbers of Type 2 AGN 
at such high X-ray luminosities are small so the results could be 
affected by small number statistics. 

In the highest X-ray luminosity bin, the AGN are expected 
to have a rel atively narrow Ed dington ratio distribution based on 
the models of lAird et alj ( l2013h . These AGN therefore likely corre¬ 
spond to high-mass black-holes hosted in massive galaxies. If AGN 
were preferentially found in star-forming galaxies on the main se¬ 
quence, these massive galaxies would be expected to have higher 
star formation rates. The high X-ray luminosity bins may therefore 
be dominated by a larger fraction of quiescent hosts for the Type 
1 AGN, which would naturally explain the lower 850^m fluxes. 
X-ray luminous Type 1 AGN could also represent a skewed pop¬ 
ulation of highly accreting black-holes seen during a brief phase 
of black-hole growth. As such, their host galaxies may be under- 
massive relative to their black-hole masses, and this could once 
again explain the lower dust luminosities and therefore star forma¬ 
tion rates in these sources. 

The Type 2 AGN show a more marked increase in average 
submillimeter flux with X-ray luminosity with the most X-ray lu¬ 
minous AGN having higher star formation rates. As can be seen in 
Fig.H for a flux-limited sample, the median redshift increases with 
increasing X-ray luminosity. However, we have shown in Section 
I4.2l that the 850/im flux does not evolve significantly with redshift 
in the Type 2 AGN population. Therefore the trend observed in 
Fig.[6]is unlikely to be driven by redshift. Evidence of higher star 















































Table 2. Stacked 850/^m fluxes for the Type 1 and Type 2 AGN in three bins in hard X-ray luminosity. Detections above 3cr are shown in bold and the number 
of sources in each bin is given in brackets. 



Type 1 AGN 

Type 2 AGN 

43 .5<logio(L2-io/ erg s“^)<44.0 
44 .0<logio(L2-io/ erg s“l)<44.4 
logio(L 2 -io/ erg s“l)>44.4 

0.7±0.3 (N=102) 
1.0±0.2 (N=168) 
0.7±0.2 (N=146) 

0.5±0.2(N=114) 
0.6±0.2 (N=94) 
1.3±0.3 (N=50) 


formation rates i n more X-ray luminous AGN was also found by 
iLutz et aP ( l2010h who argue for two modes of star formation in the 
population - secular and merger-driven. The most X-ray luminous 
Type 2 AGN could therefore represent a transition to a merger- 
driven scenario where star formation and black hole accretion are 
more tightly coupled. Fig. also suggests that there may be sys¬ 
tematic differences in the host galaxy properties of the most X-ray 
luminous Type 1 and Type 2 AGN with the luminous Type 2 AGN 
occupying more highly star-forming galaxies. Alternatively, it is 
possible that the Type 2 AGN have higher 850^m fluxes simply as 
a result of a higher dust mass relative to their Type 1 counterparts 
of comparable luminosity. 

Overall, it is difficult to establish a direct causal connection 
between AGN activity and star formation when studying average 
properties of AGN host galaxies over a wide range in luminosity, 
redshift and host galaxy type. The situation is further complicated 
by the fact that the timescales for star formation and AGN activity 
are very diff erent so AGN variability can wipe out any un derlying 
correlations jHickox et al.l20l4lGabor & Boumaudll2013tl . Indeed 
the AGN that are individually detected at >3.5 (t at 850/^m (Table 
|T]l have a range of hard X-ray luminosities and would populate all 
three X-ray luminosity bins considered in this analysis. However, 
at least in the Type 2 AGN population there appears to be evidence 
for a transition to more highly star-forming hosts and/or higher dust 
mass at the highest X-ray luminosities. 


4.4 Obscuration 

AGN classifications into Type Is and Type 2s as described in Sec¬ 
tion are one way to potentially separate unobscured and ob¬ 
scured AGN. However, other obscuration measures have also pre¬ 
viously been used in the literature, many of these tracing obscura¬ 
tion on potentially very different physical scales. A key aim of this 
work is to understand if any of these obscuration measures are cor¬ 
related with the 850/im flux and therefore star formation rates of 
AGN. Having established trends in the 850/im flux with both red- 
shift and X-ray luminosity in the preceding sections, we will aim 
to disentangle these effects from the effects of obscuration on the 
cold dust emission. 


4.4.1 Type 1 versus Type 2 AGN 

There are 428 z > 1 Type 1 AGN and 211 z > 1 Type 2 
AGN in our sample as classified based on their spectroscopic 
properties/hest-fit SEDs (see Section [IS. The median redshifts 
are 1.78 and 1.45 and the median hard X-ray luminosities are 
logio(L2-io)=44.2 and 44.0 respectively. In Appendix C we also 
show both the X-ray luminosity and redshift distributions for the 
Type 1 and Type 2 populations. The Type 2 population clearly 
shows a paucity of sources at the highest redshifts. However in 
Section l4^ we have already established that redshift trends are un¬ 
likely to be significant within the population. The stacked 850/im 
flux densities for the Type Is and Type 2s are 0.8±0.1 mjy and 


0.6±0.1 mJy respectively. In Section l43l we have seen that the 
high luminosity Type 2 AGN have higher submillimeter fluxes but 
at lower luminosities, the Type 1 AGN 850/im fluxes are slightly 
higher. Given that the X-ray luminosity distribution for the Type 
2 AGN sample peaks at lower luminosities than the Type 1 AGN, 
this helps explain the small difference in the overall stacked fluxes 
observed here. 


4.4.2 X-ray Hardness Ratio 


The X-ray hardness ratio (HR) is used to characterise the colour 
index associated with the X-ray spectrum of an AGN and there¬ 
fore the level of absorption seen in the X-ray. The HR is de¬ 
fined as [H + S)/ [H — S) where H represents the counts in 
the XMM-Newton hard X-ray band and S represents the counts 
in the XMM- Newton soft X-ray band. An HR value of —0.2 has 
been used by iBrusa et^ 1 I 2 OIOI) to separate obscured and unob¬ 
scured AGN and following this, we split our AGN sample into 
two HR bins. Only those AGN that are detected in both the soft 
and hard X-ray bands are used for this analysis. Our sample com¬ 
prises 129 HR^ —0.2 obscured AGN (34 Type 1 and 95 Type 2) 
and 346 HR< —0.2 unobscured AGN (263 Type 1 and 83 Type 
2). The average 850/im fluxes are very similar for the obscured 
and unobscured samples as shown in Table - 0.7±0.3 mJy at 
HR^ —0.2 and 0.7±0.1 mJy at HR< —0.2. In Appendix C we 
also show the redshift and luminosity distributions for the two sam¬ 
ples, which are very similar. Various previous studies have consid¬ 
ered the dependence of cool dust emission from AGN on X-ray 
obscuration estim ates parametrized by the neutral hydrogen col- 
umn density (e.g. |S tp^enset_alJ|200j; jshaoet ^l20ld : lLutz et all 


l20ld : [Rosario et ^ 20121 : Merloni et alj 2014 ). In most of these 

studies, no significant correlation is seen between cool dust emis¬ 
sion and the dust column towards the central accreting black-hole. 
Thus the dust responsible for the X-ray absorption is probably not 
directly associated with star formation processes. Indeed most X- 
ray studies place warm absorbers on much smaller spatial scales 
than typical star-forming regions in galaxies. 


4.4.3 Optical-to-infrared colours 

Optical to infrared flux ratios have been proposed as a way of iden- 
tifyi ng luminous obscu red quasars that are not in X-ray surveys 
(e.g. lHickox et alj2007h and recent results suggest that far infrared 
fluxes and therefore star formation r ates could be corr elated with 
such colour measures of obscuration l lChen et alj|2015l) . We there¬ 
fore split our AGN sample into two (r-4.5/im) colo ur bins using 
the thr esho ld of (r-4.5/rm)=6 .1 (Vega) used in both iHickox et all 
( l2007h and IChen et al.l ( l2015h . The red sample has 303 AGN (65 
Type 1 and 238 Type 2) at a median redshift of 1.56, whereas the 
blue sample has 391 AGN (359 Type 1 and 32 Type 2) at a me¬ 
dian redshift of 1.74. The stacked 850/im fluxes are presented in 
Table[3]and are very similar for the two populations - 0.7±0.1 mJy. 
Our result, that red and blue AGN have very similar submillimeter 





































Table 3. Summary of stacked 850/rm fluxes for subsamples of AGN selected according to various obscuration parameters. Detections above 3cr are shown in 
bold. 


Parameter Investigated 

Sample 

N 

Z 

logio(L 2 io/erg s ^) 

Ssso (mJy) 

Spectroscopic/Photometric Classification 

Type 1 

428 

1.78 

44.2 

0.8±0.1 


Type 2 

271 

1.45 

44.0 

0.6±0.1 

Hardness Ratio 

Type 1 & 2 HR > -0.2 

129 

1.60 

44.2 

0.7±0.3 


Type 1 & 2 HR < -0.2 

346 

1.61 

44.3 

0.7±0.1 

Optical/Infrared Colour 

Type 1 & 2, (r-4.5)>6.1 (Vega) 

303 

1.56 

44.1 

0.7±0.1 


Type 1 & 2, (r-4.5)<6.1 (Vega) 

391 

1.74 

44.2 

0.7±0.1 


Type 1 & 2, (r-4.5)<6.1 (Vega) 
loglo(L 2 -lo) >44.4 

124 

2.15 

44.7 

0.7±0.2 


Type 1 & 2, (r-4.5)>6.1 (Vega) 
loglo(L 2 -lo) >44.4 

60 

2.31 

44.7 

1.2±0.3 


Type 1 & 2, (r-K)>5 (Vega), logio(X/0)>l 

149 

1.59 

44.2 

0.7±0.2 


Type 1 & 2, (r-K)>5 (Vega), logio(X/0)>l 
logio(L 2 -io) >44.4 

31 

1.96 

44.7 

1.6±0.4 

SED Fit Extinction 

Type 1 & 2 E(B - V)agn > 0.3 

348 

1.50 

44.1 

0.7±0.1 


Type 1 & 2 E(B - V)agn < 0.3 

351 

1.80 

44.2 

0.7±0.1 


Type 1 & 2 E{B - V)gal > 0.3 

298 

1.63 

44.2 

0.9±0.1 


Type 1 & 2 E(B - V)gal < 0.3 

366 

1.60 

44.1 

0.7±0.1 


properties, is apparently at odds with that of IChen et al.l ll201.'5h . 
However, these authors have investigated differences in red and 
blue AGN at the highest quasar luminosities. iLutz et al.l l l201(]h also 
find that trends with obscuration start to become more significant at 
high X-ray luminosities and as illustrated in Fig. differences in 
Type 1 and Type 2 AGN submillimeter fluxes only become signifi¬ 
cant at high luminosities. We therefore only select those AGN with 
logio(L 2 -io) >44.4 and again split these into two samples based 
on the (r — 4.5/rm) colours. The red sample has 60 AGN (18 Type 
1 and 42 Type 2) and the blue sample has 124 AGN (116 Type 1 
and 8 Type 2). The 850/rm fluxes in the two bins are once again 
given in Table [3] An excess in the 850^m flux is now observed in 
the red AGN sample although this is only significant at the ~1.3 (t 
level given the sample size. 

Another interesting population of red, dusty AGN ha ve also 
been identified in the COSMOS field bv iBrusa et al.l l l2010tl . These 
authors use the (r — K) colour and X-ray to optical flux ratio 
to select obscured AGN. Follow-up observations of this popula¬ 
tion demonstrate that t hey are hosted in starbur st / main-sequence 
star-forming galaxies l lBongiomo et ai] l2014tl with strong evi¬ 
dence for powerful outfl ows affecting the ionised gas ( iBrusa et al.l 
l2015l:ICresci etal.l2015h . consistent with these dusty quasars being 
caught during a brief evolutionary phase when both star formation 
and black hole accretion are synchronously occurring. The most 
X-ray luminous red AGN in the COSMOS held, XID2028, was 
independently selected in a wide-held search for re ddened, high- 
redsh ift quasars in the UKIDSS Large Area Survey jBanerii et al.l 
l2012h . The population of hyper-luminous re ddened quasars also 
shows evidence for very high dust lumin osities l Banerh^ alj2014h 
and large reservoirs of molecular gas jFeruglio et alj|2014lf . sug- 
gesting that the dust obscuration is, at least partially associated 
with star formation in the quasar host. This population of mod¬ 
erately obscured or reddened AGN appears to domi nate the AGN 
luminosity function ev en at very high luminosities dBanerii et al.l 
l2015l : ll^cv et al.ll2015h . These are therefore candidates where we 
might expect excess submillimeter emission relative to the rest of 
the AGN population. 

The brightest red AGN X1D2028 lies in a region of poor 


sensitivity in our 850^m wide map, with crsso ~ 3.4 mjy and 
is therefore not presen t in the 850/rm cat alogue. XID18, another 
obscured quasar from iBrusa et^ l l2015h . is associated with an 
850/im bright source in Table [T] We once again employ a stack¬ 
ing analysis in order to test whether such red obscured quasars are 
associated with high cold dust luminos ities and therefore high star 
formation rates. Specifically, we follow iBrusa et al.l ( 1201 Oh and iso¬ 
late a population of AGN with (r-K)>5 (Vega) and logio(X/0)>l, 
where (X/O) represents the hard X-ray to optical r-band flux ra- 
tio and the monoc hromatic r-band flux is estimated as described in 
IBrusa et ^ ( l2010h . Our sample comprises 149 AGN (14 Type 1 and 
135 Type 2) at a median redshift of 1.59. The average 850^m flux 
is 0.7±0.2 mJy and there is no strong evidence for an excess in the 
submillimeter fluxes compared to the various other obscured popu¬ 
lations that we have been investigating above. However, once again 
isolating only the X-ray luminous sub-set of this red population 
(logio(L 2 -io) >44.4), we now find that the stacked 850^m flux of 
1.6 ± 0.4 mJy to be more than 2x higher at a significance of ~2cr, 
once again suggesting that an increase in cold dust emission with 
obscuration only becomes significant at high luminosities. While 
we have put forward the interpretation that the cold dust emission 
is directly tracing star formation in these high luminosity, obscured 
AGN samples, it is also possible that the more highly obscured X- 
ray luminous AGN simply have higher dust mass, which is respon¬ 
sible for driving the observed 850/rm emission trends. However, as 
dust mass in galaxies is closely tied to their gas content, large dust 
masses generally imply large gas reservoirs and therefore presum¬ 
ably also higher star formation rates. 


4.4.4 E(B- V) from SED Fitting 

In iBongiomo et alj l l2012ti . extinctions have been derived for the 
COSMOS AGN sample from SED fitting to the extensive multi 
wavelength photometry. Separate extinction values were fit for 
the AGN and galaxy components. As expected the Type 2 AGN 
are dominated by high extinctions towards the AGN component 
whereas the Type 1 sources are fit by relatively modest extinc- 














































tions to the AGN component. We begin by separating the 2 > 
1 AGN sample based on these AGN component extinction val¬ 
ues. Two samples are created with E{B — y)AGN < 0.3 and 
E{B — y)AGN > 0.3 respectively. The low extinction sample 
comprises 349 Type 1 AGN and only 2 Type 2 AGN and has a 
median redshift of 1.80 and a mean E{B — 1^)agn = 0.04. The 
high extinction sample has 269 Type 2 AGN and 79 Type 1 AGN 
with a median redshift of 1.50 and a mean E{B — = 1 . 3 . 

The average 850/rm fluxes are 0.7±0.1 mjy in both the low and 
high extinction bins respectively and there is therefore no evidence 
for a dependence of star formation rate on these AGN component 
extinction values. 

We now consider the extinction values fit to the galaxy com¬ 
ponent instead, using a threshold of E(B — 17)gal=0.3 to sepa¬ 
rate the AGN into two samples. The reliability of SED-based ex¬ 
tinction estimates for AGN host galaxies has not been explored in 
detail thus far and should clearly be interpreted with caution, par¬ 
ticularly for Type 1 AGN. Before using these extinction values we 
therefore look for systematic biases in the two populations. Firstly, 
we c onfirm using the re st-frame monochromatic UV-luminosities 
from ILusso et^ ( 1201 dl for the Type 1 AGN, that there are no 
systematic differences in the UV-luminosity between the two ex¬ 
tinction bins (logio(Luv/ergs“^cm“^Hz“^)=29.7 in both bins). 
Next we calculate the rest-fr ame 5.8/jm lumino sity directly from 
the infrared photometry from iBrusa et al.l (1201 Oh and look for dif¬ 
ferences in the mid infrared luminosity between the two extinction 
bins. A systematic difference in the rest-frame 5.8/rm luminosity 
between host galaxy extinction bins would indicate that AGN split 
by their host galaxy extinction have different mid infra-red SEDs, 
which could reflect differences in covering factor and viewing an¬ 
gle. Once again, we find no notable difference in the rest-frame 
5.8/rm luminosity as a function of extinction, with average values 
of logiofLs.s/erg s“^)=45.1 in both bins. 

Having confirmed that the Type 1 AGN split by host extinction 
do not show systematic differences in either their UV or IR lumi¬ 
nosities, we now proceed to look at the average 850/im flux in these 
two bins. The average fluxes are quoted in Table[3and are 0.7±0.1 
mJy and 0.9±0.1 mJy in the low and high extinction bins respec¬ 
tively. The mean E(B-V)gal values in these bins are 0.06 and 0.6 
respectively. The excess 850/im flux in the high extinction bin is 
only significant at the 1 . 4(7 level. The median stellar masses and 
hard X-ray luminosities of the AGN in these two extinction bins 
are very similar - logio(M*/M 0 )=lO. 8 O and 10.85 respectively and 
logio(L2-io/ergs“^)=44.14 and 44.16 respectively. The luminos¬ 
ity and redshift distributions of the various E(B-V) selected samples 
are also shown in Appendix C and demonstrate that the samples 
overlap significantly in both X-ray luminosity and redshift space. 

Overall, we find that trends in the 850^m flux with differ¬ 
ent measures of obscuration are weak when the AGN sample is 
considered as a whole. However, we do find that these trends be¬ 
come more significant at the highest luminosities where a tran¬ 
sition from secular to merger-driven co-evolution may be occur¬ 
ring. Much larger samples over wider fields are required to confirm 
these results and investigate the overlap between obscured popula¬ 
tions selected using different classification methods. Starting with 
mid infrared rather than X-ray selected parent AGN samples would 
also help in getting better statistics on luminous, obscured AGN 
where high levels of obscuration cause the AGN to drop out of 
flux-limited X-ray samples. These studies are now becoming pos¬ 
sible wi th very large mid in frared surveys such as the WISE All-Sky 
Survey dWright et al.ll2010l) . 
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Figure 7. Reds hift versus ste l lai' ma ss of '^114,000 UltraVISTA galaxies 
at 2 : > 1 from [MuzzineLalJ i2Q13h compared to the Type 1 and Type 2 
AGN jB ongioni^^nr20 1 2h . In the case of the UltraVISTA galaxies, the 
coloured contours represent the density of points on the plot with red con- 
tou rs conesponding to t he highest density. UltraVISTA galaxies that satisfy 
the IPonlev et alJ d2012h IRAC AGN selection criteria are also shown and 
occupy the full redshift-mass plane. On the other hand, the X-ray AGN host 
galaxies correspond to the most massive end of the UltraVISTA galaxy stel¬ 
lar mass distribution at 2 > 1. 


4.5 AGN Versus Non-AGN 


Finally, we would like to understand if there is evidence for cool 
dust emission being preferentially associated with AGN activity or 
whether all galaxies of similar mass and at similar redshift, show 
similar dust emission properties. Stellar masses are available for 
both the Type 1 a nd Type 2 AGN samp les from spectral energy 
distribution fitting dBongiorno et al.ll2012h . The stellar masses as a 
function of redshift are shown in Fig. [T] As described in Section 
1431 we also have a control sample of galaxies from the UltraV- 
ISTA survey with p hotometric redshifts and stellar masses from 
iMuzzin et alj ( l2013h . The stellar mass versus redshift distribution 
for ^114,000 galaxies selected from this sample is also shown in 
FigJU, In both cases, the stell ar masses have been computed using 
the iBruzual & Chariot ( l2003h stellar population synthesis models 
with a Chabrieil ( 1200; ) initial mass function and assuming expo¬ 
nentially declining star formation histories so the two stellar mass 
estimates shoul d be comparable. Match ing the UltraVISTA galaxy 
catalogue to the lBongiorno et aT] ( l2012h catalogue, we find that for 
the Type 2 AGN sample, where the emission is dominated by the 
galaxy SED and the contribution from the AGN template is negligi¬ 
ble, the stellar masses in the two catalogues agree to within 0.1 dex 
with a median difference of 0.08 dex in the masses. In the analy¬ 
sis that follows, we separate sources into stellar mass bins of width 
0.5 dex so for our purp oses , differences in stella r masses between 
iBongiomo et alJ ( 1201 2 l) and iMuzzin et al.l ( l2013h are negligible. 

From Fig. |7] it is apparent that the AGN host galaxy masses 
are at the top end of the distribution of masses in the UltraV¬ 
ISTA galaxy population. We therefore begin by selecting both 
samples of galaxies and AGN to be in the stellar mass range 
10.0<logio(M*/MQ) <11.5 and require the AGN stellar masses 
to have flag 0 in the catalogue of iBongiorno et all ( |2012|) . which 
















































indicates that these stellar mass estimates are robust. In the case of 
the Type 1 AGN, these ^ag 0 sources correspond to those where the 
galaxy contribution at 1 is >10% so the mass of the galaxy can 
effect ively be constrained jMerloni et al.l l2010l : iBoneiomo et al.l 
l2012h . We also restrict our analysis to only those AGN with errors 
of <0.1 dex on their stellar masses. In order to remove any con¬ 
tamination from AGN in the non-AGN galaxy sample, we cross¬ 
match the U ltra VISTA galaxies to t he XMM-Newton point-source 
catalogue of ICanpelluti et all (l2009ll and remove any galaxies as¬ 
sociated with an X-ray point source within a matching radius of 
8”. Even within the mass range selected for our AGN versus non- 
AGN galaxy comparison study, the mass distribution for the AGN 
is more strongly skewed towards higher masses than the corre¬ 
sponding galaxy distribution. As it is well known that more massive 
galaxies have higher star formation rates, this can introduce biases 
into our analysis. Starting with the non-AGN, we therefore match 
both the mass and redshift distributions of the galaxies to that of the 
AGN by randomly sampling from the galaxies in (Ilog(M*) bins of 
width 0.1 and 6z bins of width 0.1. The number of AGN to non- 
AGN in each mass bin is fixed to be the same and we therefore 
ensure that the median stellar mass of the non-AGN galaxy sample 
is the same as that of the AGN sample - logio{Mt/M q) =10.9 in 
both cases. The median redshifts are 1.66 for the AGN and 1.68 for 
the non-AGN. A potential concern is that highly obscured AGN 
that are not detected in the X-ray, still remain in the non-AGN 
galaxy sa mple. However , visua l inspection of the best-fit galaxy 
SEDs from iMuzzin et al] ( l2013h confirms that the majority of these 
galaxies do not have power-law type AGN dominated SEDs in the 
Spitzer IRAC bands. From our sample of 782 mass-matched non- 
AGN galaxies, we fi nd that <5% satisfy the IRAC AGN selection 
criteria presented bv IPonlev et alj (l2012h. As illu s trated in Fig.|3 
Ultra VISTA galaxies that satisfy the lDonlev et alj j2012l) AGN se¬ 
lection criteria, occupy the entire redshift-mass plane and do not 
dominate at high masses. In this high-mass regime, the Ultra VISTA 
galaxies are therefore unlikely to suffer from heavy contamination 
from highly obscured AGN. 

With the mass and redshift distributions of the inactive galax¬ 
ies and AGN now matched, we proceed to calculating their average 
850/im fluxes. The non-AGN galaxy sample has a higher 850/im 
flux compared to the AGN - 1.01±0.09 mjy versus 0.72±0.12 
mjy. The median and mean 850/rm fluxes for both the AGN and 
non-AGN are also consistent with each other confirming that the 
stacked properties are not being dominated by a few outliers. To 
investigate this excess in more detail, we now split our galaxy and 
AGN sample into three stellar mass bins and consider the Type 1 
and Type 2 AGN populations independently. The average submil¬ 
limeter fluxes in the three stellar mass bins are presented in Table 
|4] and illustrated in Fig. [8] While the Type 2 AGN 850/rm fluxes 
increase as a function of stellar mass, the Type 1 AGN fluxes re¬ 
main fairly constant. The redshift and luminosity distributions of 
the AGN in all three stellar mass bins can be seen in Appendix 
C. In the case of both the Type I and Type 2 AGN, the highest 
stellar mass bins peak at higher X-ray luminosities consistent with 
more massive galaxies having more powerful accreting black holes. 
Combined with the results in Fig.|2 this suggests that the trends in 
submillimeter emission with mass are equivalent to the trends with 
X-ray luminosity, with Type 1 AGN showing little evolution and 
high luminosity/high-mass Type 2 AGN being more submillime¬ 
ter bright. Interestingly however, the galaxy stacked fluxes are in 
general comparable to or higher than those seen in the AGN popu¬ 
lation. For Type 2 AGN this appears to be the case over pretty much 
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Figure 8. Average 850/rm flux versus stellar mass for both Type 1 and Type 
2 AGN and compai'ed to a sample of non-AGN galaxies selected from the 
UltraVISTA sui-vey. 


the entire mass range probed. Non-AGN hosts therefore appear to 
be more highly star-forming than the X-ray AGN. 

These results are somewhat at odds with similar analyses con¬ 
ducted using Herschel data where evidence for excess star for¬ 
mation is seen in AGN host galaxies relative t o non-AGN host 
galaxies of similar mass (e.g. lSantini et al.l|2012h . The excess dis¬ 
appears when considering star-forming galaxies only leading the 
authors to conclude that AGN are preferentially hosted in star¬ 
forming galaxies with the same gas supply likely fuelling both pro- 
cesses. However , we n ote that in contrast to these investigations, 
iBongiomo et alj ( 12012h find that AGN are preferentially hosted in 
redder, more quiescent galaxies. IBongiomo et alj ( l2Q12h have al¬ 
ready pointed out that their SED-fitting method, by correcting for 
dust extinction to both the AGN and host galaxy components and 
disentangling the UV-luminous AGN emission from that from star 
formation, leads to lower sta r form ation rate estimates for AGN 
hosts relative to Santini et al.l ( l2012h . The AGN template used by 
ISantini et alj ll2012h is a Sevfert template. While this is clearly ap¬ 
propriate for low-luminosity AGN, for the high-luminosity XMM- 
COSMOS AGN investigated here, it may not account for all the 
UV emission from the AGN, therefore leading to an over-estimate 
of the star formation rate. 

The higher average 850^m flux of the non-AGN galaxy 
sample is primarily driven by the most massive galaxies with 
logic (M*/Mq) >11.0 as can be seen in Fig. [8] These most mas¬ 
sive galaxies have extremely red SEDs with very little flux at rest- 
frame UV wavelengths. Many of these ultra-red galaxies would not 
be selected in flux-limited catalogues generated in bluer bands such 
as the r or i-band and could be one rea son for the different conclu- 
sions from our analysis and that of e.g.lSantini_et_^ ( l2012h . Inter¬ 


preting the higher 850^m fluxes of the iMuzzin et al j ( l2013h mas¬ 
sive galaxies as signifying that they have higher star formation rates 
relative to AGN, requires us to assume that the average dust proper¬ 
ties are similar in both populations. Inspection of th e best-fit non- 
AGN galaxy SEDs at logic(Mt/M©) >11.0 from iMuzzin et al.l 
( l2013h however confirms that these galaxies are almost always fit 
by highly star-forming galaxy templates with significant extinction 
(median Ay = 0.9), lending credence to th e idea that a s ignifi- 
cant fraction of the massive galaxies found bv iMuzzin et alj ( l2013ll 



























































Table 4. Stacked 850/im fluxes for the Type 1 and Type 2 AGN and a sample of non-AGN galaxies selected from UltraVlSTA in three stellar mass bins. 
Detections above 3cr are shown in bold and the number of sources in each bin is given in brackets. 



Type 1 AGN 

Type 2 AGN 

Non-AGN UltraVlSTA Galaxies 

!0.0<logio(M,/MQ)<10.5 

lO.5<logio(M*/M0)<ll.O 

ll.O<logio(M*/M0)<11.5 

0.6±0.3 mJy (N=35) 

1.0±0.3 mJy (N=97) 
1.0±0.2 mJy (N=91) 

0.3±0.4 mJy (N=22) 
0.4±0.2 mJy (N=91) 
0.9±0.2 mJy (N=69) 

0.6±0.2 mJy (N=98) 
0.8±0.1 mJy (N=329) 
1.4±0.1 mJy (N=355) 


at these redshifts may indeed correspond to starburst galaxies with 
very high star formation rates and therefore submillimeter fluxes. 

We have shown that AGN host galaxies appear to have lower 
850^m fluxes compared to a mass-matched sample of non-AGN 
galaxies. It is important to highlight that photometric redshifts as 
well as stellar masses for both the AGN and non-AGN samples are 
clearly model dependent and affected by the choice of templates 
used in the SED fitting. Stellar mass estimates are also likely to 
be less reliable for Type 1 AGN where the host galaxy is less visi¬ 
ble. However, as can be seen in Fig.[^ it is the difference between 
the submillimeter fluxes of Type 2 AGN hosts and non-AGN that 
is more significant. Different model choices in SED fitting could 
perhaps partially explain the different conclusions reached regard¬ 
ing the nature of AGN host galaxies in the literature. However, an¬ 
other key issue is that many previous studies have selected non- 
AGN galaxy samples from optical survey data, which are less sensi¬ 
tive to star-forming galaxies with significant dust extinctions. High- 
mass, high-redshift, red populations of non-AGN galaxies that are 
matched in mass to the most luminous AGN are now selectable us¬ 
ing new deep infrared surveys such as UltraVlSTA that cover rea¬ 
sonably wide fields. Due to their rarity, these massive red galaxies 
are only present in wide-field extragalactic surveys such as COS¬ 
MOS and due to the red colours of this population, they are often 
not present in optical flux-limited surveys. These galaxies are also 
the most difficult to obtain spectra of. Understanding their true na¬ 
ture - i.e. whether they correspond to massive dusty starbursts - is 
now vital in order to be able to relate them to the AGN host galaxies 
which have equivalently high stellar masses. 


5 CONCLUSIONS 

We have studied the 850/rm properties of X-ray selected AGN in 
the COSMOS field using new data from the SCUBA-2 Cosmol¬ 
ogy Legacy Survey. The 850/rm data covering ~2 deg^ presents an 
opportunity to study the cold dust emission from a statistically sig¬ 
nificant number of luminous AGN as a function of properties such 
as luminosity, stellar mass and obscuration. We begin by searching 
for AGN that are detected at >3.5cr in our 850^m data in a “high 
sensitivity” (ctsso <2 mJy) region covering 0.89 deg^. There are 
360 submillimeter sources over this area and we identify that 19 
of these are X-ray luminous AGN. These submillimeter galaxies 
have S85 o= 4-10 mJy. The 19 submillimeter detected X-ray AGN 
comprise 6 spectroscopically confirmed broad-line AGN, 1 spec¬ 
troscopically confirmed narrow-line AGN, 3 photometric redshift 
AGN that are best fit by Type 1 AGN templates and 9 photomet¬ 
ric redshift AGN that are best fit by Type 2 or galaxy templates. In 
terms of X-ray luminosity, these AGN fill the luminosity gap be¬ 
tween the very X-ray luminous quasars that ha ve been individually 
followed up at submillimeter wavelengths (e.g. lStevens et alj2005h 
and sub-mm detect ed X-ray AGN from much smaller area, deeper 
X-ray surveys (e.g. [Alexander et aI]|2005l : IWang et al.ll2013l) . 

We construct inverse-variance weighted stacked 850^m im¬ 
ages of different sub-samples of z > 1 AGN in the X-ray lu¬ 


minosity range 43<logio(Lx)<45, utilising a hard X-ray (2-10 
keV) flux-limited population of 699 radio quiet AGN (428 Type 
1 and 271 Type 2). The average 850/rm fluxes are measured from 
a 100 such bootstrap images and for the entire AGN population, 
S85 o= 0.71±0.08 mJy. The 850/rm flux directly traces the dust lu¬ 
minosity at 2 ; > 1 due to the effects of the negative k-correction so 
the 850^m fluxes can be averaged over relatively wide redshift bins 
at these high redshifts. Under the assumption of an isothermal SED, 
the 850/rm flux can therefore be used as a proxy for star forma¬ 
tion rate. Using independent estimates of the star formation rate in 
the COSMOS AGN sample estimated via SED fitting to the optical 
through infrared photometry, we demonstrate that AGN with higher 
star formation rates do, on average, have higher 850/im fluxes. In 
particular, AGN with SFR> lOMQyr”^ have 850/rm fluxes that are 
'~3.3x higher than those in AGN with SFR<10MQyr“^. We then 
study the dependence of this 850/rm flux on various AGN proper¬ 
ties. In particular, we reach the following conclusions: 

• The 850/rm fluxes of both Type 1 and Type 2 AGN show little 
evolution with redshift and the Type 1 and Type 2 AGN have av¬ 
erage submillimeter fluxes that are consistent with each other over 
the entire redshift range probed in this study. 

• We study the dependence of 850/rm emission and there¬ 
fore star formation on AGN X-ray luminosity. We find that the 
850/rm fluxes are relatively constant for z > 1 Type 1 AGN at 
logj^Q(L 2 -io) > 43.5. However, in the case of the z > 1 Type 2 
AGN, we find that the 850^m flux increases with X-ray luminosity. 
This may suggest a transition from secular to merger-driven evolu¬ 
tion in the Type 2 AGN population at the highest X-ray luminosities 
(logio(L 2 -io) > 44.4) probed in this work, which could be re¬ 
sponsible for the tighter coupling seen between cold dust emission 
(tracing star formation) and X-ray luminosity (tracing black hole 
accretion) at these luminosities. Our results also indicate that high- 
luminosity Type 1 and Type 2 AGN have systematic differences 
in their observed 850/rm flux densities, which could reflect sys¬ 
tematic differences in their host galaxy properties with the Type 2 
AGN occupying more highly star-forming galaxies. Alternatively, 
the higher 850/im flux densities in high luminosity Type 2 AGN 
could simply arise as a result of these AGN having higher dust 
mass relative to Type 1 AGN of comparable luminosity. 

• Having explored trends with redshift and luminosity, we study 
the dependence of submillimeter flux on obscuration, relying on 
various different obscuration measures. There is no significant de¬ 
pendence of the 850/rm flux on X-ray hardness ratio or the (r- 
4.5/im) colour when averaging over all X-ray luminosities. When 
selecting the most X-ray luminous AGN (log]^g(L 2 -io) > 44.4), 
we find that the red AGN have submillimeter fluxes that are ~1.7x 
higher than the blue AGN, although given the sample size, the re¬ 
sult is only sig nificant at the ~1_3 (t level. 

• Following iBrusa et^ ( l2010h . we select a population of red 
(r — K) AGN with high X-ray to optical flux ratios. Once again 
these AGN have very similar submillimeter properties to the rest 
of the AGN population when averaging over all X-ray luminosi¬ 
ties. However in the highest X-ray luminosity bin (log]^g(L 2 -io) > 













44.4), these red AGN have submillimeter fluxes that are more than 
a factor of 2 higher than the average. Despite the small samples 
in this single field, the results are still significant at the ~2(j level 
and consistent with a more marked dependence of submillimeter 
emission on obscuration for the most X-ray luminous AGN. 

• Finally, the z > 1 AGN are compared to a A'-band selected 
population of galaxies showing no evidence for AGN activity and 
that are carefully matched in stellar mass and redshift to the AGN 
sample. On account of the Ff-band selection, this non-AGN galaxy 
sample includes populations of extremely red galaxies at the high- 
mass end, many of which have very little flux at rest-frame UV 
wavelengths and would not be selected in optical flux-limited sur¬ 
veys. Type 2 AGN show an increase in submillimeter flux with stel¬ 
lar mass which is equivalent to the increase seen as a function of 
X-ray luminosity. Over almost the entire mass range probed how¬ 
ever, non-AGN have higher 850/im fluxes compared to the Type 2 
AGN and in the highest stellar mass bin - 11.0 <logio < 11.5 - 
the non-AGN also have higher submillimeter fluxes compared to 
the Type 1 AGN. The best-fit SED templates for the massive non- 
AGN galaxies indicate that they are dusty galaxies with high star 
formation rates and a median Av of 0.9. Such massive, red non- 
AGN galaxies have eluded discovery in previous surveys that have 
either covered smaller areas or been flux-limited in the bluer optical 
bands. Due to their red colours, many of these galaxies also remain 
spectroscopically unconfirmed. Understanding their true nature is 
important to be able to relate them to their AGN counterparts of 
equivalent mass. 

We have demonstrated the unique sensitivity and resolution of 
the new SCUBA-2 850^m data provided by the SCUBA-2 Cosmol¬ 
ogy Legacy Survey, which has allowed us to probe rest-frame wave¬ 
lengths corresponding to the Rayleigh Jeans tail in the SEDs of 
luminous X-ray AGN in the high redshift Universe. At these wave¬ 
lengths the 850/rm emission provides a unique probe of the cold 
dust luminosity and dust mass in AGN samples. Although the sam¬ 
ple size used in this work has limited the significance of some of our 
results, nevertheless we find some evidence for a stronger coupling 
between AGN luminosity, dust luminosity and obscuration at high 
X-ray luminosities, where we may be seeing a transition from a 
secular to a merger-driven regime in galaxy evolution. Larger sam¬ 
ples of AGN selected at multiple wavelengths are clearly required 
to establish more significant trends in the star formation properties 
of AGN host galaxies across a wide range in luminosity, redshift 
and mass. This will become possible when similar analyses are ex¬ 
tended to cover all the fields imaged by the S2CLS. 
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APPENDIX A: X-RAY LUMINOSITIES & ABSORPTION 
CORRECTIONS 

When deriving hard X-ray luminosities for our AGN sample in Sec¬ 
tion we have neglected the effects of absorption. Significant 
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Figure Al. Hard X-ray luminosities deriv ed in this paper (Sectionl2.2t ver¬ 
sus those derived using Chandra data bv iBiightman et aU l2014l) and in¬ 
cluding the effects of absorption. Filled circles correspond to AGN at 2: > 1 
which are used for all the stacking analysis considered in this work. No sys¬ 
tematic bias is seen in our X-r ay luminosities compar ed to the absorption 
corrected X-ray luminosities of iBrightman et all l2014ll . confirming that ab¬ 
sorption has a negligible impact on the derived luminosities for the majority 
of our sample. 


absorption in a large fraction of our sample would systematically 
bias our X-ray lu minosity estimates. As illustrated in figure 5 of 
lAird et all ( 1201 5h however, hard X-ray luminosities are only likely 
to be significantly biased for absorbing columns of Nh > 10^® 
cm~^. A subset of the COSMOS AGN considered in this work 
also have deep X-ray observations obtained using the Chandra ob¬ 
servatory and where the data is of sufficient quality to enable de- 
tailed spectral fi t ting i ncluding different absorption components. 
iBrightman et alj ( l2014h have conducted exactly such a study and 
derived absorption corrected hard X-ray luminosities as well as col¬ 
umn densities for th e X-ray AGN. We there fore match our sample 
of X-ray AGN to the lBrightman et al] ( l2014ll sample. There are 217 
AGN in common between the two catalogues and where the red- 
shifts in both catalogues also agree. In Fig. lAll we compare the hard 
X-ray luminosities derived in this paper ( without absorption correc - 
tion) to the hard X-ray luminosities from iBrightman et alj (l2014ll . 
As can be seen, most of the AGN lie on the x = y relation and 
th e median difference bet ween our X-ray luminosities and those 
of iBrightman et all ( l2014h is 0.06 dex. When considering the de¬ 
pendence of submillimeter flux on hard X-ray luminosity, we stack 
sub-samples in bins of 0.4-0.5 dex in X-ray luminosity. Hence we 
conclude that any biases in our X-ray luminosities due to a lack of 
absorption correction, are unlikely to significantly affect the trends 
observed in this paper. 


APPENDIX B: PROPERTIES OE 850/im DETECTED 
X-RAY AGN 

In this section we present the properties of our sample of 19 X- 
ray AGN that are associated with 850/im sources in the COSMOS 
field. Table lB ll nresents the stellar masses and various obscuration 
parameters investigated in this work for each of the 19 AGN. Ta¬ 
ble |B^ presents photometry from Spitzer MIPS (24/im), Herschel 
PACS (100/im & 160/im) and SPIRE (250, 350 and 500/im). In 


























Table B3. Best-fit SED parameters for individually detected AGN with at 
least two detections in the Herschel bands. 


XID 

omir 

Tdust/K 

logio(LFiR/L0) 

logic (Mduat/M©) 

13 

1.21±0.07 

15±2 

11.41±0.16 

9.08±2.41 

139 

1.55±0.09 

22±1 

11.85±0.09 

8.56±0.34 

160 

1.97±0.06 

37±2 

12.70±0.05 

8.58±0.07 

270 

1.84±0.04 

31±2 

12.57±0.07 

8.73±0.07 

278 

1.54±0.08 

32±5 

12.34±0.14 

8.52±0.08 

353 

1.63±0.07 

20±2 

11.69±0.15 

8.70±0.15 

53922 

1.59±0.04 

21±2 

11.80±0.10 

8.62±0.21 

54440 

1.79±0.07 

27±4 

12.35±0.14 

8.79±0.10 

60070 

1.54±0.10 

25±4 

12.17±0.18 

8.79±0.11 


Fig.im we show the best-fit SEDs for the AGN that are detected 
in at least two of the Herschel bands. We fit a power-law of the 
form S\ oc A“ to model the mid infrared emission from th e AGN 
at A > 3/im, together with a single temperature greybody dCasevI 
l2012n with a fixed dust emissivity index of /3 = 2. The best-fit mid 
infrared power-law slope, a, dust temperature, far infrared lumi¬ 
nosity (60-300/im) and dust mass are given in Table lB3l In Fig. lB2l 
we also show the correlation between the far infrared luminosity 
and the 250 and 850/im fluxes for these individually detected AGN 
as well as the correlation between the 250 and 850^m fluxes and 
the dust mass. 


APPENDIX C: REDSHIFT & X-RAY LUMINOSITY 
DISTRIBUTIONS 

Here we present the redshift and hard X-ray luminosity distribu¬ 
tions for the various sub-samples that have been stacked in this 
analysis. These distributions are shown in Fig. m an d ica 








Table Bl. Summary of the properties of the 19 X-ray AGN that are associated with 850^m sources. Note: the hard X-ray luminosities and redshifts of these 
sources are presented in Tabled 


XID 

logio(M*/M0) 

HR 

(r-4.5^m) Vega 

E(B-V)agn 

E(B-V)gal 

13 

11.26 

-0.58 

5.44 

0.1 

0.0 

18 

11.17 

-0.21 

7.68 

0.6 

0.0 

139 

10.08 

-0.37 

6.61 

0.6 

0.2 

160 

10.61 

-0.09 

5.18 

0.8 

0.2 

246 

9.12 

0.42 

8.59 

1.3 

0.3 

250 

11.11 

- 

5.02 

0.1 

0.5 

270 

11.74 

- 

8.13 

0.3 

0.3 

278 

11.26 

- 

9.13 

2.6 

0.5 

353 

11.24 

-0.14 

6.28 

0.0 

0.4 

402 

10.54 

-0.24 

6.57 

0.3 

0.4 

415 

9.92 

- 

6.18 

0.3 

0.5 

469 

10.44 

- 

6.00 

0.3 

0.1 

10675 

11.16 

- 

5.37 

1.0 

0.1 

10809 

10.09 

- 

6.56 

0.4 

0.0 

30182 

9.62 

-0.6 

3.82 

0.3 

0.0 

53922 

9.25 

-0.16 

6.66 

2.3 

0.2 

54440 

11.18 

- 

9.07 

0.7 

0.5 

60070 

11.32 

- 

7.99 

0.4 

0.5 

60490 

10.66 

- 

6.57 

0.3 

0.0 


Table B2. Summaiy of Spitzer MIPS and Herschel fluxes of the 19 X-ray AGN that are associated with 850^m sources. Note: the SCUBA-2 850/im flux 
densities of these sources are presented in Tabled 


XID 

S24/mJy 

Sioo/mJy 

Si6o/mJy 

S25o/mJy 

Ssso/mJy 

Ssoo/mJy 

13 

0.91 ±0.06 

- 

- 

13.8±2.8 

34.4±5.2 

23.9±4.8 

18 

0.51±0.02 

- 

- 

<9.2 

<10.6 

<12.2 

139 

2.27±0.37 

18.4±2.9 

55.7±3.0 

59.6±6.5 

45.7±5.9 

31.3±5.7 

160 

1.59±0.12 

29.3±2.1 

75.2±3.2 

69.3±7.0 

51.6±6.2 

29.0±5.4 

246 

0.88±0.02 

- 

- 

<9.0 

<12.1 

<12.6 

250 

0.28±0.02 

4.9±1.2 

- 

<9.2 

<10.6 

<12.2 

270 

0.87±0.02 

13.7±1.6 

31.3±3.1 

53.3±6.0 

48.8±6.0 

22.5±5.0 

278 

0.56±0.02 

- 

14.9±3.5 

18.3±3.2 

21.9±4.1 

16.2±3.6 

353 

0.40±0.02 

- 

12.5±2.9 

16.8±3.1 

<16.1 

<16.3 

402 

0.28±0.09 

- 

13.3±3.3 

<14.6 

<15.9 

<13.6 

415 

0.13±0.02 

- 

- 

< 9.2 

<10.6 

<12.2 

10675 

0.09±0.02 

- 

- 

<9.2 

<10.3 

<12.6 

53922 

1.24±0.03 

10.7±1.8 

35.2±3.7 

36.4±4.9 

34.7±5.2 

- 

54440 

0.39±0.02 

- 

15.7±3.8 

27.1 ±4.1 

27.2±4.7 

<16.3 

60070 
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Figure Bl. Best-fit SEDs for the 9 AGN in Table lB2l with at least two detections in the Herschel bands. 















Figure B2. Top: Correlation between the far infrared luminosity estimated from full SED fitting (proxy for star formation) and the 250 and 850/i.m fluxes of 
individually detected AGN. Bottom: Correlation between the dust mass estimated from full SED fitting and the 250 and 850/im fluxes of individually detected 
AGN 
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Figure Cl. Redshift (left) and X-ray luminosity (right) distributions for various sub-samples of AGN compared in this work. The top panels show AGN 
senarateH as Tvne I and Tvne 7 the second nanels show AGN senarated based on hardness ratio the third nanels show AGN senarated on the basis of ontical- 
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Figure C2. Redshift (left) and X-ray luminosity (right) distributions for the Type 1 and Type 2 AGN split into the thi'ee stellai' mass bins considered in Section 

EH 








































































